In [39]:
import sys, os
import numpy as np
import matplotlib.pyplot as plt

# Import the aims module
from soma import aims
# the brainplot package
import colorado as cld

print(sys.version)
3.6.9 (default, Jan 26 2021, 15:33:00) 
[GCC 8.4.0]

Anatomist snapshot

In [40]:
%%capture 
# The previous magic caputres this cell output (namely Anatomist logs)

# Import anatomist and set the environment to use it here
import anatomist.api as anatomist
os.environ["QT_API"] = "pyqt5"
%gui qt

# create an Anatomist session and capture input (so it does not appear here)
a = anatomist.Anatomist();
In [41]:
# load data (the SliceableObject)
t1mri = a.loadObject("data/subject01.nii")

# create an Axial window in anatomist
w = a.createWindow("Axial", geometry=[1200, 350, 500, 500])
t1mri.addInWindows(w)

# get a snapshot from anatomist
%matplotlib inline
img = cld.anatomist_snatpshot(w)
plt.imshow(img);
In [42]:
## View meshes
lwhite = a.loadObject("data/subject01_Lwhite.mesh")
rwhite = a.loadObject("data/subject01_Rwhite.mesh")
w3d = a.createWindow("3D")
w3d.addObjects([lwhite, rwhite])

plt.imshow(cld.anatomist_snatpshot(w3d))
Out[42]:
<matplotlib.image.AxesImage at 0x7f9ddeaf8470>

Meshes

In [43]:
meshR = aims.read('data/subject01_Rhemi.mesh')
meshL = aims.read('data/subject01_Lhemi.mesh')
cld.draw([meshL, meshR])

Buckets

In [44]:
# load a bucket with aims
bck_map=aims.read('data/roi_Bucket.bck')
a = bck_map[0]
cld.draw(a)
In [45]:
# load a bucket from a bucket txt file (needs sulci_isomaps)
import sulci_isomap as siso

dx, bucket, _ = siso.files_utils.parse_bucket_file('data/LAbby_New.bck')
cld.draw_numpy_buckets([bucket])

Volumes

In [46]:
vol = aims.read('data/subject01.nii')
cld.draw(vol, max_points=10000, th_min=950, th_max=1000)
In [ ]: